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We investigate the tunability of electrostatic coupling between solid state quantum dots as building 
blocks for quantum bits. Specifically, our analysis is based upon two-dimensional electron systems 
(2DEG) and depletion by top gates. We are interested in whether the Coulomb interaction between 
qubits can be tuned by electrical means using screening effects. The systems under investigation are 
analyzed numerically solving the Poisson equation in 3D via relaxation techniques with optimized 
algorithms for an extended set of boundary conditions. These include an open outer boundary, 
simulation of 2DEG systems and dielectric boundaries like the surface of a physical sample. The 
results show that for currently lithographically available feature sizes, the Coulomb interaction 
between the quantum bits is weak in general due to efficient screening in the planar geometry of 
2DEG and top gates. The evaluated values are on the order of 1 fieV. Moreover, while it is not 
possible to turn off the qubit interaction completely, an effective tunability on the order of 50% is 
clearly realizable while maintaining an intact quantum bit structure. 

PACS numbers: 03.67.Lx, 85.35. Be, 85.35.Gv 



I. INTRODUCTION 
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Quantum dot structures based upon localized regions 
of charge play an important role in investigating quan- 
tum mechanical effects at the nanoscale. 1 ' 2 ' 3 Despite the 
short dccohercnce times, considerable progress is made 
in utilizing the charge degree of freedom for coherent 
quantum processes with possible application in quantum 
information. 4 ' 5,6 In particular, electrostatically-defined 
quantum dots based on depletion of two-dimensional 
electron systems (2DEG), play an important role due to 
their variability as the charge "drop" in the dot can be 
controlled in shape and size by applying well-controlled 
gate voltages. 

In this paper we consider a linear array of double dots 
as shown in Fig. 1. With different arrangement of voltage 
gates, quantum dots are formed and arranged in pairs to 
define charge qubits. Every pair is defined with a tun- 
nel contact (quantum point contact, QPC) between the 
two dots, which can be tuned by an external gate volt- 
age. With the QPC in a weakly transparent regime, the 
charges are well localized in their respective quantum dot, 
and the dot occupation is characterized by a well-defined 
number of electrons. The system can then be tuned to 
the energetically degenerate case, where an extra electron 
can reside on either quantum dot symmetrically The 
dynamics of such a system is well-described by the two 
qubit states \L) and \R), with the extra electron either 
on the L(eft) or on the R(ight) dot. These two states 
form the natural basis to define the two charge qubit 
configurations easily controlled by the gates that define 
the quantum dots and the QPC' 

As different dot pairs are used to define charge qubits 
via local gate voltages, their states are well defined 
and controlled experimentally. However, it is crucial to 
be able to control the interactions between neighboring 
qubits if one is to implement the quantum gate oper- 



FIG. 1: Linear array of double quantum dots. Charge is 
exchanged between the dots in each pair but different pairs are 
only coupled electrostatically. Defining gates are not shown. 



ations required in quantum computation. We are thus 
interested here in the energetics of different charge config- 
urations in the double-dot array and the tunability of the 
coupling between charge qubits via purely electrostatic 
means. Notice that including the interaction with neigh- 
boring qubits results in an effective dipole-dipolc interac- 
tion (the common ZZ interaction in solid state qubits in 
a spin notation). This causes the states in the qubit to 
align according to their neighbors, and misalignments to 
be energetically costly. We investigate the energetics of 
different qubit configurations and how one can stabilize 
different states by purely electrostatic means. 

For typical nanoscopic devices with many (or at least 
a few) electrons in each of the electrostatically-defined 
regions, the charge distribution and the major energy 
scales are described to a good approximation by classical 
electrostatics. 9 Due to the strong electric fields generated 
by segregating charge in a 2DEG, the Coulomb energy is 
the dominant energy scale. Thus, it is important to know 
the electrostatics of the system if one expects a good 
quantitative description thereof. We present here a de- 
tailed study of the electrostatic coupling in the quantum 
dot array of Fig. 1. We find that one can substantially 
modify the qubit couplings by proper gate geometry and 
voltage control (~ 50%). However, we also find that it is 
not possible to completely turn off the coupling among 
dots, as one would need in order to allow for total decou- 
pling of the system. Our studies utilize an efficient three- 
dimensional electrostatic code that uses higher-order grid 
relaxation and fast Fourier transform algorithms. This 
allows for accurate and realistic modelling of the quan- 
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turn dot array of interest here, as well as a variety of other 
structures. Effects such as depletion in a 2DEG due to 
shallow etch 10 or deep etch 11 have been investigated with 
our code and the results agree well with experimental 
data and other simulations built on Schrodinger Poisson 
solvers. 12 

As we mentioned above, we are interested in the tun- 
ability of the coupling between charge qubits. Electro- 
statically tunability is directly related to changes in the 
geometry of the system, i.e. the regions were there is 
charge, since a constant geometry of electrically isolated 
regions implies a constant capacitance matrix and con- 
stant electrostatic coupling. Application of gate voltages 
produces changes in the lateral charge depletion on the 
2DEG system, drastically changing the electric field line 
distribution around the area of interest. In the quan- 
tum dot structures analyzed here, this includes field lines 
that emerge out of the dots through the surface, as well 
as throughout the semiconducting material including the 
depleted 2DEG region. Notice, however, that depiction 
of charge in the 2DEG already suggests that the local 
electric field is strong with dominant field lines vertical 
to the 2DEG layer. Small local changes of charge or rear- 
rangements are likely not to affect the interdot coupling 
strongly. This conclusion is not trivially reached, due to 
the complexity of the dot array and gate geometry which 
includes subtle screening effects. Still, this intuitive pic- 
ture turns out to be correct, as will be demonstrated in 
detail in the following. 



II. SIMULATION OF 2DEG 

2DEG systems are "volatile" by design, i.e. the electric 
charge can be depleted or accumulated by external means 
using appropriately fabricated conducting leads acting as 
voltage gates. Thus the 2DEG regions containing elec- 
tric charge density change dynamically with the voltage 
patterns applied to the gates. For a flexible and reliable 
numerical analysis it is thus necessary that the 2DEG 
boundary is treated dynamically. The plane of the 2DEG 
may contain several independent regions of charge that 
can be addressed individually. Moreover, well isolated 
regions of charge have a floating potential, and this sit- 
uation remains even if weak tunnel junctions are present 
because of charge quantization. These considerations in- 
troduce a set of specific boundary conditions related to 
2DEG systems. They are readily realized in our numer- 
ical code and explained in detail in [13]. Moreover, as 
from a numerical point of view we are dealing with a fi- 
nite system, the potential on the outer boundary floats 
freely and the system must be solved self-consistcntly. 
For more detail in that respect the reader is also referred 
to the Appendix. 

The primary boundary condition for the 2DEG as it 
turns out is surprisingly simple local update prescription 
built in easily with the relaxation sweeps: 13 when coming 
across a 2DEG grid point, say point i, with the 2DEG 



FIG. 2: (Electric potential V(z) in a parallel plate geometry 
to demonstrate depletion in 2DEG - the material below the 
surface (to the right of the ai layer) is considered uniform 
with a dielectric constant of e. The charge densities on the 
three planes shown are ai, a p (p for the positive charge of 
the donor layer) and 02, respectively. Vo is the total potential 
difference across the whole stack related to the pinning of the 
Fermi level at the surface. 



layer kept from the exterior at the potential V2DEG i then 
proceed as follows - first, update the potential V] locally 
as in homogenous space. Second, if V" cw > V2DEG- take 
Vi = V2DEG otherwise accept the calculated V/"™. The 
result is that the regions where V} = V2DEG contain 
electronic charge while the regions in the 2DEG with 
Vi < V2DEG ar< 2 depleted. The basic idea underlying 
this prescription is simply that for the electronic charges 
qi < in the 2DEG it is energetically favorable to go 
to a potential that is higher then V^deg- The effect is 
that charges gather in these regions until the potential 
equilibrates at V^deg- 

The gross features of the charge distribution in a real- 
istic sample with a uniform 2DEG system buried in close 
vicinity to the surface can be well understood by a par- 
allel plate capacitor arrangement as shown in Fig. 2. o\ 
is the charge density at the free surface. A distance di 
below lies the donor layer ( "delta doping" ) with a doping 
density a v = enp, followed by the 2DEG separated by a 
distance d 2 from the donors. The variables d\ and d 2 are 
parameters that can be adjusted to account for the finite 
width of each of the three layers. The electrical voltage 
Vo of the surface with respect to the 2DEG is fixed as 
a consequence of the pinning of the Fermi energy at the 
surface and makes up for the difference in chemical po- 
tential due to surface states. For the case of Ga(Al)As, 
for example, the Fermi level is pinned to mid-gap, 1 and 
so the potential Vo is chosen to be fixed at Vo ~ —0.7 V 
(note that the 2DEG is considered grounded). 

The charges at the surface and in the 2DEG rearrange 
and equilibrate when the temperature is high enough, 
typically at room temperature. This leads to the ener- 
getically most favorable charge distribution that shares 
the charges available from the donor layer between the 
surface and the 2DEG such, that the electrical field above 
the surface and below the 2DEG of the sample is zero. 
The remaining adjustable parameter in a real system is 
then the distance d\ which can be reduced by etching. 
The critical distance d\ for depletion of the 2DEG layer 
below is given by df lt = —Voeeo/eno which is indepen- 
dent of the distance d 2 (note that Vo is negative here). 

With clearly separated layers of uniform charge, this 
capacitor model can be expected to provide a good de- 
scription of charge distribution. The quantum mechani- 
cal effects are implicitly present in the effective parame- 
ters d\ and d 2 as well as in effects such as the pinning of 
the Fermi energy at the surface. To this extent, this semi- 
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classical description of the system provides an accurate 
picture. 9 



III. ENERGETICS OF DOUBLE DOT QUBITS 

Double dot systems are created in the 2DEG via deple- 
tion through negatively biased metallic top gates on the 
surface of the structure, e.g. fabricated by lithographic 
means. The typical system analyzed here is shown to- 
gether with its numerical results in Fig. 3. The relatively 
large number of gates defining the structure keeps the 
system flexible. However, the exact number of gates and 
their respective arrangement are arbitrary to a certain 
extent and may be simplified by fewer gates in a specific 
experimental setup. 

Every double dot system (qubit) in our design has al- 
together six gates defining the dots, including a plunger 
gate for every dot, and two gates that tune the tunnel 
barrier (QPC) between the two dots defining the qubit. 
One of the two QPC gates for a pair of dots may eventu- 
ally be merged with one of the gates confining the dot. 

The set of top gates chosen can be seen in panel (a) 
of Fig. 3. The charge distribution shown with biased 
gates clearly outlines their geometry. The red shading 
(negative charge density) results from the negative bias 
voltages applied. Thus the more negative the bias, the 
larger the accumulated negative charge on that gate. 

The potentials on the top gates were chosen as follows 
depending on the local degree of depletion that was re- 
quired: —2.5 V on the plunger gates for every dot, —5.2 V 
for the QPC between a pair of dots, —2.4 V for the tran- 
sition region between pairs of quantum dots, except for 
the center transition whose gate potential was alterna- 
tively also set to OV. For the other confining potentials 
a value of —4.5 V was taken. The remaining (geomet- 
rical) parameters are: pinning of the surface potential 
Vq = —0.75 V, dielectric constant e = 12, grid spacing 
h = 12 nm, di = di = 36 nm, and no = 4 • 10 16 m" 2 , re- 
sulting in a charge distribution ratio between the surface 
and the 2DEG of about 2:1. 

In the setup of Fig. 3, the quantum dots contain about 
200 electrons. In the following, the two adjacent qubits 
in the center of Fig. 3 are considered in more detail, and 
the outer qubits allow us to examine the effect of qubit 
pair interaction. The plunger gate for individual dots was 
fine tuned self-consistently at the end of the initialization 
process to produce half-integer occupancy on the dot. 
With this it is then an easy numerical task to achieve 
equivalent integer charge distributions for a single qubit 
with one extra electron symmetrically on either of the 
dots forming the qubit labelled as the two qubit states \L) 
and \R). More specifically, with the gate voltages frozen 
to their initial configuration for half-integer occupation of 
the dots, the charge configuration can now be altered to 
the possible different integer value charge configurations, 
shifting half of an electron in a quantum dot to the other 
or vice versa. 



FIG. 3: (Color online) Simulation of 2DEG with top gates and 
applied gate voltages. The center gates (labeled V c in panel 
a) extend in between the qubit pair such that the 2DEG un- 
derneath the gate V c can be depleted at will (see panels b 
and c). The color coding for charge distributions is consis- 
tently taken such that red corresponds to negative and blue 
to positive charge density (in gray scale all shades are to be 
considered red except for the dark gray regions in panel d). 
a) Outline of top gates by showing the charge distribution 
on these gates. The 1 /im bar shows chosen length scale, b) 
Charge distribution in the 2DEG with the center gate at the 
same voltage as the gates separating the other qubits. c) Same 
as (b) but with the voltage on the center gate lifted such that 
additional charge accumulation is allowed in the 2DEG un- 
derneath the gate, d) Differences in the charge distribution 
for the setup in (c) resulting from small variations on the 
plunger gate voltages chosen such that the total change of 
charge within one dot is ± e/2. e) Free electrostatic energy of 
the different charge configurations shown in (d) taken relative 
to the first configuration #1. 

As indicated in Fig. 3, the gate in the center is intro- 
duced for the purpose of screening the qubit-qubit in- 
teraction through depiction or accumulation of charge in 
between the qubits, thus introducing the change of ge- 
ometry of the charged regions necessary for tunability of 
the qubit interaction. The two resulting configurations 
analyzed are shown in panels (b+c) for the gate voltages 
V c = —2.4 V and V c = V, respectively. For the setup 
in panel (b), the different qubit states resulting from sin- 
gle electron hopping on the two qubits in the center are 
presented in the panels (d). For better visual contrast, 
the difference in the charge configuration Aq (r) in the 
layer of the 2DEG with respect to the initial half-integer 
setup is shown. When summing up Aq (r) over a sin- 
gle dot, then the total difference in charge AQi for dot 
i is |AQj| = e/2 up to numerical precision. Note that 
the influence on the outer qubits, i.e. the qudot pair on 
the very left and the very right as seen in panels (b+c) 
is minimal, and no appreciable charge rearrangement oc- 
curs there due to the charge alterations in the two middle 
qubits. 

From the numerical evaluation of the electrostatics, we 
obtain the potential of the quantum dot array together 
with the total charge on every dot. For a calculated set of 
inital values q a and V a for the charge and the potential 
on the individual dots, the electrostatic free energy is 

evaluated using AW a b = J V (q' ; V g ) dq' . The vector V g 

la 

stands for a set of external potentials held at constant 
voltage. For the case of weak changes in the systems 
electrical geometry, the total change in free energy for 
different charge configurations is given by 13 

AW ab -&q\ {V a + V b ) , (1) 
where Aq = q b — q a . For a constant geometry (capaci- 
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tance matrix) this equation holds exactly. However, in 
the relevant cases here where gate voltages do change 
the geometry, namely the shape of the quantum dots as 
well as the remaining extended regions of charge in the 
2DEG, Eq. (1) only holds to first order in small varia- 
tions of the geometry Yet this is definitely the case for 
the single-electron hopping processes considered here. 

The qubit-qubit interaction energy for the configura- 
tions in panel (d) is then evaluated by means of Eq. (1) 
with the results shown in panel (e) of Fig. 3. The quan- 
tum dots have been intentionally designed not to be ex- 
actly the same, resulting in small relative charge varia- 
tions (maximum deviation of two electrons between the 
dots). This variation is also reflected in the slight as- 
symetry seen in panel (e) between the ideally equivalent 
cases shown in the panels (d) . For the setup in panel (b) 
with the charge depleted under the center gate, the qubit- 
qubit interaction energy is given by AE ~ 1.17/xeV. By 
lifting the potential on the center gate from V c — —2.4 V 
to V c = OV, thus allowing extra charge to screen the 
qubit-qubit interaction (see middle of panel c), the in- 
teraction energy is reduced to AE ~ 0.66 fj,eV which is 
about half the former value. 

The quantum dot spacing is about 300 nm from center 
to center of charge within a qubit, and 430 nm between 
neighboring qubits. Considering the initial half- integer 
configuration as neutral, i.e. the negatively charged 
qudots as screened by the environment, then a simple 
classical estimate of the energy scale for the qubit-qubit 
interaction by shifting half of an electron charge in be- 
tween the qubit gives 65 /ieV, which is two orders of mag- 
nitude larger than the actual realization of dots formed 
by depletion. This strong suppression of the interqubit 
interaction is clearly due to the large metallic top gates 
which efficiently screen the slightly rearranged charge 
configurations closeby. 

A similar calculation was performed with the center 
gate stretching only half-way through the qubit-qubit in- 
termediate regions. The results arc shown in Fig. 4. The 
average distance between the qudots is slightly different 
from the previous case (320 nm intraqubit vs. 430 nm in- 
terqubit separation) . The confining potentials on the top 
gates are slightly weaker, allowing larger dots with about 
300 electrons each. The remaining parameters are the 
same as in the previous case. 

The dots are completely symmetric in this case, re- 
flected in the energy of the equivalent charge configura- 
tions in panel (d) plotted in panel (e). Again, the qubit- 
qubit interaction energy for each configuration in panel 
(d) was evaluated using Eq. (1), and results shown in 
panel 3(e). For the setup in panel (b) with the charge 
depleted all along the center gate, the qubit-qubit in- 
teraction energy is given by AE ~ 0.71 fxeV. By rais- 
ing the potential on the center gate from V c = —2.5 V 
to V c = —0.75 V, the interaction energy is reduced to 
AE ~ 0.29 /ieV which is less than half the former value, 
but still clearly present. The simple classical estimate 
similar to the previous example with the extra charge in 



the qudots positioned around the center of charge leads 
to 71 /ieV, again two orders of magnitude larger than the 
energies in the actual realization. 



IV. CONCLUSIONS 

The tunability of qubit interaction based on electro- 
static gate screening has been analyzed numerically. Our 
program uses an efficient relaxation algorithm, and ef- 
ficiently describes general 3D structures, including the 
ability to handle an extended set of boundary conditions, 
specifically appropriate for the treatment of 2DEG sys- 
tems and including dielectric effects. 

The typical systems analyzed are based on 2DEG ge- 
ometries derived from GaAs/GaAlAs structures and de- 
pletion via top gates. For a chain of qudot pairs - the 
qubits - the region in the 2DEG between qubits is de- 
pleted at will, allowing for tunable screening and control 
of the nearest neighbor qubit interaction. The tunability 
derived for realistic gate voltage changes is on the order 
of 50%. The interaction energy can thus be clearly mod- 
ulated. However, we find that it cannot be turned off 
for moderate voltages, as it would be desirable for the 
on/off-switching of the interaction of ideal qubits. We 
believe that this general trend does not only hold for the 
few systems analyzed in this paper. The tunability based 
on charge rearrangement and screening of Coulomb in- 
teraction on charged-localized quantum dots is definitely 
possible, but limited in efficiency It may indeed be nec- 
essary to implement dynamic pulse control techniques, 
such as the "bang-bang" approach or others discussed 
for qubits recently. 15 
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FIG. 4: (Color online) Same as Fig. 3 but with the center 
gate (labeled V c in panel a) extending only halfway through 
from top to bottom. This structure has also slightly increased 
quantum dot spacing and size from that in Fig. 3 . The color 
coding for charge distributions is the same, with red for neg- 
ative and blue for positive charge density (in gray scale all 
shades are to be considered red except for the dark gray re- 
gions in panel d). See description of panels (a-e) in caption 
of Fig. 3. 
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APPENDIX - RELAXATION ON THE GRID 
AND OPEN BOUNDARIES 

The algorithm used in this paper is based on a higher- 
order relaxation algorithm for a finite physical system 
as described by the authors in an earlier paper. 16 Ma- 
jor modifications are significantly optimized update of 
the outer boundary, and an extended set of additional 
boundary conditions mainly related to 2DEG systems. 
As later ones represent dielectric media it is also clear 
that a proper treatment of surfaces requires a proper in- 
corporation of the dielectric constant. A detailed descrip- 
tion can be found in reference [13]. 

Solving the electrostatics in 3D for a finite system 
with an open outer boundary may not allow for a sim- 
ple boundary condition such as constant potential on the 
outer surface. With uniform grid spacing, one is bound 
to a finite total grid size, which in 3D is even more strin- 
gent. Therefore the region under investigation cannot 
assume constant potential on the outer boundary. How- 
ever, by having direct access to the local charge distribu- 
tion throughout the grid during relaxation, one can calcu- 
late the potential on the outer boundary. The naive way 
of summing up charge over distance Vi ~ J2j Qi/dij, how- 
ever, is computationally expensive, to the extent that the 
actual relaxation process becomes much faster compared 
than the self-consistent update of the outer boundary. 

The reason for this lies in the folding of the charge ver- 
sus the relative distance. The structure of the sum sug- 
gests the use of fast Fourier transform (FFT) algorithms 
in 3D 1 ' where the folding is replaced by a single prod- 



uct for each Fourier component. The fact that FFT au- 
tomatically implies periodic boundary conditions which 
may introduce artificial effects can be circumvented by 
doubling the system size in every dimension. 18 Despite 
the largely enhanced number of grid points in 3D, the 
resulting update of the outer boundary is drastically two 
orders of magnitude faster than the original way of car- 
rying out the sum in real space. 

Typical materials such as Ga(Al)As or AI2O3 have a 
large dielectric constant of order e = 12. This strongly 
affects the screening of charges such as those in the 2DEG 
or surface charges, and it is absolutely crucial to include 
the dielectric properties in these calculations. The nu- 
merical implementation goes along the lines of algorithms 
found in the literature, 19 ' 20 and it is discussed in detail 
in [13]. The main idea is that the local relaxation of the 
potential must be done under the constraint that there 
is no free charge in bulk dielectric media. This leads 
to the incorporation of the dielectric constant into the 
weights of the potential average over neighboring sites, 
ynew = Y^ jeNN tjQj / Y, je KN e j ■ 0ne should be aware 
that the dielectric constant €j has to be evaluated mid- 
way between the grid point i under consideration and the 
neighboring point j. Finally, note that the subsequent 
calculation of the source qi from the potential without 
weighting it by the dielectric constant gives the induced 
charge at dielectric boundaries (note that in uniform di- 
electric media, the dielectric constant drops out in the 
averaging process for the local update of V z new and so 
the induced charge equals zero by construction). 
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